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Abstract. Neglecting the second order corrections in weak lensing measurements can lead to 
a few percent uncertainties on cosmic shears, and becomes more important for cluster lensing 
mass reconstructions. Existing methods which claim to measure the reduced shears are not 
necessarily accurate to the second order when a point spread function (PSF) is present. We 
show that the method of Zhang (2008) exactly measures the reduced shears at the second 
order level in the presence of PSF. A simple theorem is provided for further confirming our 
calculation, and for judging the accuracy of any shear measurement method at the second 
order based on its properties at the first order. The method of Zhang (2008) is well defined 
mathematically. It does not require assumptions on the morphologies of galaxies and the 
PSF. To reach a sub-percent level accuracy, the CCD pixel size is required to be not larger 
than 1/3 of the Full Width at Half Maximum (FWHM) of the PSF, regardless of whether the 
PSF has a power-law or exponential profile at large distances. Using a large ensemble (> 10 7 ) 
of mock galaxies of unrestricted morphologies, we study the shear recovery accuracy under 
different noise conditions. We find that contaminations to the shear signals from the noise 
of background photons can be removed in a well defined way because they are not correlated 
with the source shapes. The residual shear measurement errors due to background noise are 
consistent with zero at the sub-percent level even when the amplitude of such noise reaches 
about 1/10 of the source flux within the half-light radius of the source. This limit can in 
principle be extended further with a larger galaxy ensemble in our simulations. On the other 
hand, the source Poisson noise remains to be a cause of systematic errors. For a sub-percent 
level accuracy, our method requires the amplitude of the source Poisson noise to be less than 
1/80 ~ 1/100 of the source flux within the half-light radius of the source, corresponding to 
collecting roughly 10 source photons. 
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1 Introduction 

Weak gravitational lensing has been widely used as a direct probe of the mass distribution 
of our Universe on different scales, including large scale structure, clusters, galaxies, etc. [1]. 
Not only is the physics of lensing well understood in the context of General Relativity, but the 
lensing effect can also be straightforwardly measured using the shapes of background galaxy 
images [2-4]. Currently, one of the main challenges in this field is about how to accurately 
recover the cosmic shear field from galaxy shapes [5-8]. This is difficult due to the large 
galaxy shape noise, the involvement of the point spread function (PSF), the presence of the 
photon noise, the pixelation effect, etc.. There have been many literatures focusing on this 
particular topic [9-25]. 

In all of the practical shear measurement methods proposed so far, there is a common 
assumption: the cosmic shear is small, therefore the second or higher order terms in shear 
can be neglected. This is true for the shear field of our Universe on large scales, which is 
typically of order a few percent. However, future weak lensing survey may require shear 
measurement accuracy to be controlled below a 0.1% level [26, 27]. On arc minute angular 
scales, the second order terms can cause a systematic error of order 10% to the cosmic shear 
power spectrum [28, 29]. More importantly, the shear field by a foreground cluster can easily 
be of order ten percent or more. Neglecting second order terms in shear measurements can 
lead to significant errors on the implied cluster masses [30-37]. 

The main purpose of this paper is to further develop the shear measurement method 
of [24](Z08 hereafter) by including the second order terms in shear /convergence in the for- 
malism. This is done straightforwardly in §2, in which we show that to the second order 
in accuracy, Z08 measures exactly the reduced shears. A popular misunderstanding in the 
weak lensing community is that all of the existing shear measurement methods are already 
accurate to the second order in shear /convergence, because they all claim to measure the 
reduced shears. We show why this is not generally true in the presence of PSF. In §3, we 
provide a simple and useful theorem for judging whether shear estimators are accurate to 
the second order based on their properties at the first order. The theorem provides an easy 
way to see why the method of Z08 yields exactly the reduced shear. §4 demonstrates the 
accuracy of Z08 using a large number of computer-generated mock galaxies of unrestricted 
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morphologies in the presence of PSF and the photon noise, including both the background 
noise and the source Poisson noise. Finally, we conclude in §5. 

2 The Shear Measurement Method 

The basic idea of Z08 is to use the spatial derivatives of the galaxy surface brightness field 
to measure the cosmic shear. In a parallel paper [38], we have shown that this method is 
equivalent to measuring shears using the galaxy quadrupole moments in Fourier space, with 
an additional term correcting for the PSF. Indeed, the measurement should be carried out in 
Fourier space, in which the moments can be easily evaluated, and the PSF can be transformed 
into the desired isotropic Gaussian form through multiplications. Further more, as shown in 
the parallel paper, not only the measurement, but also the whole analytic derivation of the 
relation between the galaxy surface brightness and the cosmic shear can be worked out in 
Fourier space in a much simpler way than in real space. As will be shown below, because of 
this convenience, the relation can be made accurate even to the second order in shear. 

Let us start the discussion in real space. We define the galaxy surface brightness dis- 
tribution before lensing as fs(x ), the lensed galaxy (before being processed by the PSF) as 
fL{x L ), and the observed image as fo(x°), where x s is the coordinate in the source plane, 
and x L and aP are the positions in the image plane. We have the following relations: 



h(x L ) = fs(x S ), x s = Mx L , 

fo(x°) = J d 2 x L W(x° - x L )f L (x L ), (2.1) 



where W is the PSF, M is the lensing distortion matrix typically defined as: My = 5ij — 4>ij 
with (f>ij = 5{j — dxf /dxj being the spatial derivatives of the lensing deflection angle. 4>ij 
is often replaced by the convergence k [= (4>n + $22) /2] and the two shear components 71 
[= (0n — ^22) /2] and 72 (= $12) ■ Note that here we adopt the convention used in most other 
lensing literatures {e.g. , [11, 39]) to define the shear components and the convergence. They 
are different from the definitions in the other papers of this series [24, 25, 38] at the second 
order. Their differences are shown and discussed in the appendix. 

Z08 has shown that the following relations can be used to measure the cosmic shear: 



1 «(<9l/ ) 2 " {d2f ?) 9 )en 

= -71 ) 



2«(<9 1 /o) 2 + (d2/o) 2 + A) 9 ) en 

({difod 2 fo) g ) 
«(3i/o) 2 + (d 2 /o) 2 + A) 9 ) en 

where di denotes d/dxi, and 



-72, (2.2) 



A = ^V/ -V(V 2 /o). (2.3) 



2 

is the scale radius of the isotropic Gaussian PSF Wp defined as: 




^ W = 2^ eXP l 

(• • • ) g means taking the average over the surface brightness field of a single galaxy. (• • • ) en 
means taking the average over an ensemble of galaxies. Note that whatever the original 
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PSF is, it is always transformed into the desired isotropic Gaussian form that is defined in 
eq.(2.4). This is allowed as long as the scale radius /3 of the target PSF is somewhat larger 
than that of the original PSF. With the help of Fourier transformation, we now show how to 
improve the accuracy of eq.(2.2) to the second order in shear/convergence. Readers who are 
not interested in the mathematical derivations can jump to eq.(2.13) and eq.(2.14) for the 
main results. 

For a general 2D field /(x), let us use f(k) to denote its Fourier transformation, which 
is defined as: 



/(*) 



(2.5) 



The terms in eq.(2.2) can be written as integrations of the Fourier modes of the images 
weighted by proper functions of the wave vector as follows [38]: 



{d i fo{x)d j f {x)) g 
V/o • V(V 2 /o)' 



l fd?k 

S (2tt) 2 n%n ' 
1 fd 2 k 

k 



fo(k) 



fo(k) 



(2.6) 



g S (2vr)2 

where S is the total area of the map. Consequently, eq.(2.2) can be written in Fourier space 



as: 



(P20 — Po2) t 



where 
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Pij = j d 2 kk[ki \f (k) 



-7i, 



-72, 
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Note that D4 = P40 + 2P22 + Po&- To find out how Pij changes under lensing to the second 
order in shear /convergence, we use the following relations that are derived from eq.(2.1) and 
properties of Fourier transformation: 



d z x 



det 



ik L -(M.~ 1 S s ) f s ($S^ 



\dx S J 

|det(M- 1 )| j d 2 x s e i( - M - lkL > sS f s (x s ) 
|det(M" 1 )|/;(M- 1 fc L ), 



fo(k) = Wp(k)f L (k). 



(2.9) 



Note that in the above equations, we have assumed that the PSF is isotropic Gaussian (Wp), 
whose Fourier transformation is Wp(k) [= exp(— f3 2 \k\ 2 /2)]. The dependence of Pij on the 
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cosmic shear can be derived as follows: 



Pij = / d 2 kk\k{ W^(k)\det(m- l )\f s (Mr 1 k) 



|det(M 



-In 



d 2 k(Mk)\(Mk) J 2 Wp(Mk)f s (k) 



(2.10) 



The last step is achieved by re-defining M 1 k as k. Expanding eq.(2.10) up to the second 
order in shear /convergence, and using the fact that the intrinsic galaxy shapes are isotropic 

2 



(i.e. , ( f s (k) ) en only depends on 
show the following: 



(P20 — Po2) t 
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Note that 0( 7 3 ) refers to terms of the third or higher orders in 71,2, K. From eq.(2.11), we 
find the main results of this paper: 



(P20 — Po2) t 



7i 



2 (P20 + P02 ~ (3 2 D 4 /2) 

( P ll)en 



1 - K 

72 



(P20 + ^02 - /3 2 D 4 /2) en 1 



+ 0(7 3 ), 
+ 0(7 3 )- 



(2.13) 



The terms on the right sides of eq.(2.13) are called reduced shears. For convenience, in the 
rest of the paper, we use to represent 71,2/(1 — k) respectively. One can equivalently 
write down the above formula in real space as: 



1 (((dlfo? - (d2fo?) 9 ) 



glen 
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We have just shown that to the second order in accuracy, our shear estimators yield exactly 
the reduced shears. It is important to note that the ensemble averages in eq.(2.13) / eq.(2.14) 
should be taken before the ratios. This point is extensively discussed in [38]. 

The novelty of our result may be questioned by those who believe that the existing 
shear measurement methods already claim to measure the reduced shears. Such a judgement 
is mainly based on the fact that the lensing distortion matrix M can be re- written as: 

which shows that k changes the galaxy size, and the shape distortions are exclusively due 
to the reduced shears g\^- Therefore, the results of shape/shear measurements should only 
depend on the reduced shears, not the convergence, as long as the shape measurement is 
decoupled from the galaxy size. This is usually the case in the absence of PSF. When a 
PSF is present, however, the convergence k can modify the observed galaxy shape through 
changing the galaxy size relative to that of the PSF. Generally, for spin-2 shear estimators 1 
such as those in most other existing methods, the results can generally be written as power 
series of and k as 2 : 

5i i2 [a + bK + c{gl + gl) + ---] , (2.16) 

in which a, b, and c are coeficients that may require calibrations (a = 1 and b = in our 
method). Obviously, unless b is proved to be zero, a shear measurement method does not 
directly measure the reduced shear at the second order level. For example, if we mistakenly 
use (1/2) (P 20 - Po2)en/(P20 + Po2)en and (Pn) en / {P20 + Pm)en defined in eq.(2.11) as shear 
estimators {i.e. , neglecting the term D4 for the PSF correction), we find the following 
formulae to the second order in accuracy: 



1 {P20 ~ Pm) en _ ( h A ( p n)en n /„ , ,y. 
9 IP I P \ ~ -9l{a + 0K), —————— -g 2 {a + bK), 

2 {^20 + P02)en \^20 + P02) en 



(2.17) 



in which 



u 



1--, b=u +u >-^y^, ^ 2 kff- 

\ 2 / en \ 2 / en 

This simple example has shown that if accuracy at the second order is desired, one should pay 
attention to the possible existence of terms proportional to #1,2 k in addition to the reduced 
shears. If one plans to achieve the second order accuracy by calibrating the multiplicative 
factor a + bK using numerical simulations, it is important to realize that the factor typically 
depends on k. These troubles are not present in our method. 



3 A Useful Theorem 



Since the pre-lensing galaxy sizes are not known a priori, it would be ideal for the statistical 
expectation of a shear estimator to depend only on the reduced shears, i.e. , to be completely 
decoupled from the galaxy sizes, or the convergence k. In terms of eq.(2.16), this corresponds 
to b = 0, and the coefficients in front of other high order k terms are all zero. This is, 
however, highly nontrivial to achieve in the presence of PSF, especially when all high order 

1 According to [38], we only need to consider spin-2 shear estimators. 
2 The formula is initiated by the referee of this paper. 
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shear/convergence terms are considered. Less ambitiously, we have shown in the previous 
section that the shear estimators of Z08 indeed decouple from the galaxy sizes up to the 
second order in shear /convergence. Our calculation is quite laborious, and the situation is 
likely similar in all other shear measurement methods. Fortunately, we find a simple way of 
knowing the property of a shear estimator at the second order according to its first order 
results: 

Theorem 3.1 For spin-2 shear estimators whose statistical expectation values (in the pres- 
ence of PSF, but without noise) can be written in the form of eq.(2.16), if a = 1 (i.e. , the 
method is free of calibrations at the first order), we must have 6 = 0. 

Proof If a = 1, the statistical expectations of the shear estimators can be written as: 

gi, 2 [l + bK + c(gf + gl) + ■■■], (3.1) 

For a given set of observed galaxies, let us consider the following two cases with fixed g\ and 
g2- 

1. k = 0; 

2. and the galaxies are intrinsically larger than those in the first case by a factor 

Of (1 - Hi). 

Since the observed galaxy images in the two cases are identical, the shear measurement 
method should yield the same results. On the other hand, according to eq.(3.1), we should 
get (if accurate to the second order) 51,2 in the first case, and 51,2(1 + bn) in the second case. 
Therefore, we have 6 = 0. 

The above theorem shows that o = 1 is a sufficient (not necessary) condition for 6 = 0. 
For a being any nonzero constant, the shear estimator can be trivially rescaled to make a = 1, 
which again leads to 6 = 0. When a is a function of the galaxy morphology and size in a 
method, the first case in the above proof yields 51,20', while the second case gives 51,2(0 + 6k). 
Since a and a' can naturally have a difference of order k due to the change of the intrinsic 
galaxy sizes, the value of 6 is not necessarily equal to zero. 

The theorem can be phrased in plain words as: if a method measures exactly the shears 
at the first order without the need of calibrations, it must yield exactly the reduced shears 
at the second order level. One can immediately apply this theorem to the method of Z08 
to show why it accurately measures the reduced shear without the lengthy calculation in §2. 
The theorem is equally useful for judging the accuracy of any shear measurement method at 
the second order based on its properties at the first order. 

4 Numerical Test 

In this section, we test the accuracy of shear recovery with a large ensemble (> 10 7 ) of mock 
galaxies. In principle, eq.(2.13) / eq.(2.14) allows us to recover the cosmic shear to a sub- 
percent level accuracy. In practice, however, one needs to consider at least two important 
factors: the pixelation effect and the photon noise. The purpose of this section is to test this 
method under realistic conditions, and to show how these factors may affect the accuracy of 
our method. 

[25] has introduced useful interpolation methods to treat the pixelation effect, which 
becomes a problem for shear measurement when the CCD pixel size is comparable to the size 
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of the PSF. We find that these methods are not accurate at the sub-percent level in terms 
of shear recovery. On the other hand, in the course of this work, we surprisingly find that if 
it is only for the purpose of shear measurement, oversampling galaxy images does not seem 
necessary. At least for the PSFs used in this work, the shear estimators calculated in Fourier 
space converge quickly when the pixel size is less than about 1/3 of the FWHM of the PSF. 
This will be demonstrated in §4.2. 

The treatment of the background photon noise has also been discussed in [25]. It can in 
principle be used to deal with both the background fluctuation and its Poisson noise, since 
they are both independent of the source flux. There is no need for any modifications of the 
treatment, because the quantities that we need to measure from each galaxy are the same 
as those discussed in Z08. The treatment is indeed simple: for each galaxy, one subtracts 
from the nominators and denominators in eq.(2.13) / eq.(2.14) the contributions from the 
photon noise that are estimated from a neighboring map of pure noise [25]. Note again that 
the ratios should be taken after the ensemble averages. As will be shown in §4.3, the shear 
measurement errors due to the background noise can be removed cleanly with the treatment 
of [25]. On the other hand, we do not yet have a way to correct the shear measurement errors 
due to the source Poisson noise. In this paper, we simply perform the shear measurements 
without additional treatment of the source Poisson noise. A further development in this 
aspect will be studied in a separate work. 

4.1 General Setup 

Each of our mock galaxies is placed at the center of a 192 x 192 grid. The grid size is used 
as the length unit in the rest of this paper. The PSF has a truncated Moffat profile used in 



The FWHM of this PSF is very close to r^. We always set r c = 3rd in the simulations of this 
paper. In our shear measurement method, the PSF is always transformed into the isotropic 
Gaussian form in Fourier space before the shear measurement is carried out. The scale radius 
of the target PSF defined in eq.(2.4) is set to 0.7r d , so that the FWHM of the target 
PSF is slightly larger than that of the original PSF. is set to 12 in unit of the grid size. 
Note that the grid size is not equivalent to the pixel size. The later is always chosen to be 
an integer multiple of the grid size in this paper. 

Our mock galaxies are made of point sources that are generated by 2D random walks. 
There are at least four main purposes for this arrangement: 

1. To maximize to some extent the richness of galaxy morphologies; 

2. The lensing effect can be exactly mimicked by simply changing the positions of the 
point sources; 

3. Convolution of the galaxy image with any PSF is trivial and easy; 

4. It is extremely fast to generate such galaxies. 

Based on the above facts, we encourage everyone working in the field of shear mea- 
surement to test their methods with the random- walk-generated mock galaxies 3 . In our 

3 We are aware of the fact that our random-walk galaxies do not have certain properties of real galaxy, such 
as the sharp cusp, the long tails of the de-Vaucouleurs profile, or any systematic trends in ellipticity gradients 



the GREAT08 project (see [8] for details): 




(4.1) 
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simulations, each such mock galaxy is made of forty point sources, whose positions are deter- 
mined by forty steps of 2D random walk. Each step size is a random number between and 
1. The direction of every step is completely random in 2D. For each galaxy, the first step 
starts from the center of the grid. The ending position of the i th step is where the i th point 
source is. In the course of the random walk, if the distance of the i th point to the center 
of the grid is more than 6, we restart the i th step from the center of the map. Every point 
source of a galaxy is assumed to carry the same luminosity. The resulting post-seeing galaxy 
has roughly the same FWHM as that of the PSF. 

To test the shear recovery accuracy of our method, we use six sets of input shear values 
(7i, 72): (0.05, -0.05), (0.03, -0.03), (0.01, -0.01), (-0.01, 0.01), (-0.03, 0.03), (-0.05, 
0.05). k is fixed at 0.05. For each set of the shear values, we use about 10 7 mock galaxies 
to recover the shear. To calibrate the shear recovery accuracy quantitatively, we adopt the 
commonly used multiplicative bias m and additive bias c defined as follows: 

g measured = (1 + mjg**** + Cl , 

g rneasured = (1 + m2 )^put + ^ ^ 

For each shear component, we calculate the two bias parameters by fitting a linear relation 
between the measured and the input values of the reduced shear. 

4.2 Test the Pixelation Effect 

We set the pixel size to be a multiple of the grid size to mimic the pixelation effect. Since 
we need to use Fast Fourier Transform (FFT) in our method, the size of each galaxy postage 
stamp is chosen to be an integer power of 2 times the pixel size along both sides. The power 
integer is chosen to be the value that allows the postage stamp to cover the most area of the 
whole grid. For example, if the ratio of the pixel size to the grid size is set to 5, the postage 
stamp contains 32 x 32 pixels, because 5 x 32 < 192, and 5 x 64 > 192. Note that our choice 
of the postage stamp size is large enough for avoiding cutting off the edges of the simulated 
galaxies. The simulations in this section do not contain noise. 

Under-sampling of the galaxy images can certainly affect the shear recovery accuracy. 
We have discussed a few interpolation methods in [25] to help reduce the systematic errors 
due to the pixelation effect. The most accurate methods are found to be the so called "Log- 
Bicubic" or "Log-Spline" methods (equally good), which simply refer to performing the 
traditional "Bicubic" and "Spline" methods on the logarithm of the data instead of the data 
itself. The focus of the previous work is on the Gaussian PSF. The PSF defined in eq.(4.1) 
has a power-law decaying form when the distance to the center is large. For this type of PSF, 
it is perhaps not surprising that the Log-Bicubic method is found not to work as well as in 
the Gaussian PSF case, although it is still better than the traditional Bicubic method. On 
the other hand, fortunately, we find that shear measurement in Fourier space is quite robust 
even for marginally under- sampled images (without the need of interpolation), regardless 
of whether the PSF is a Gaussian function or a power-law form. This is demonstrated 
in fig.l and fig. 2, which shows how fast the measured value of the first shear component 



with radius. However, sharp cusp should not cause any singularity in the method because of the smoothing by 
the PSF. Neither the long tails of the de-Vaucouleurs profiles nor any systematic trends in ellipticity gradients 
with radius seems to be able to affect the accuracy of this method as long as the whole/complete galaxy image 
is captured/used. A more detailed study of these issues will be included in a future work. 
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Figure 1. This figure compares the qualities of three different ways of treating the pixclation effect. 
The input PSF is defined in cq.(4.1), with FWHMsa = 12. For each chosen pixel size L p i xe i (in 
unit of the grid size), we plot the measured first component of the reduced shear from a single galaxy 
(defined on the left side of the first half of eq.(2.13)) normalized by its measured value at L p i xe i = 1. 
The solid curve is from a direct Fourier transformation without any interpolation; the dashed curve is 
from the data interpolated by the Log-Bicubic method; the dotted line is from the data interpolated 
by the Bicubic method. In the second and last cases, the pixel size of the interpolated images is 
always set to 1. 



from a single galaxy (denned on the left side of the first half of eq.(2.13)) converges for an 
increasingly smaller pixel size. The results are shown for three methods: using the Bicubic 
interpolation method (the dotted curves); using the Log-Bicubic method (the dashed curves); 
direct Fourier transformation/shear measurement (the solid curves). Fig.l and fig. 2 are for 
the PSF defined in eq.(4.1) and the Gaussian PSF respectively. The two PSFs have the same 
FWHM. Clearly, according to the figures, the direct measurement (i.e. , no interpolation or 
any other treatment for the pixelation effect) consistently yields no more than 0.5% relative 
deviations from the correct answers as long as the pixel size is not larger than 4, which 
corresponds to about 1/3 of the FWHM of the PSF in both cases. It is useful to note that 
the ratio 1/3 between the pixel size and the FWHM of the PSF is indeed reasonable in 
practice, and used in the GREAT08 project. The Log-Bicubic method, on the other hand, 
can cause a few percent error on the shear measurement when the PSF has a power-law 
decaying form, despite its superiority in the Gaussian PSF case. 

For further checking the accuracy of direct Fourier space measurements on pixelated 
images, we use a large ensemble of mock galaxies and the input shear components given in 
the previous section to study how the values of the multiplicative and additive bias vary with 
the pixel size. Here and in the rest of the paper, we use the PSF defined in eq.(4.1). In table 
1 and 2, we show the results for six different choices of the pixel size. For each pixel size, we 
consider four choices of galaxy size, which are listed in the tables as small, medium, large, 
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Figure 2. Same as fig.l, except that the input PSF has a Gaussian form define in cq.(2.4). The 
FWHMs of the two PSFs are equal. 



largest galaxies (written as "Gal." in the tables for abbreviation). The generation of the 
medium size galaxies are described in §4.1. The small, large, largest galaxies are generated 
through the same procedures, except that the overall scales of the galaxies are multiplied by 
factors of 0.5, 2, and 4 respectively before they are convolved with the PSF. The average 
FWHMs of the post-seeing galaxies of small, medium, large, and very large sizes are equal 
to 1.02, 1.05, 1.16, and 1.43 times the FWHM of the PSF respectively. In our simulations, 
when a galaxy is generated, it is used for different shear values and pixel sizes for the purpose 
of saving time. Galaxies of different sizes are always generated using different sets of random 
seeds instead of simple rescaling. 6.5 x 10 6 mock galaxies are used to measure the biases in 
every case of the tables. 

According to table 1 and 2, for all cases with L p i xe i < 4, the multiplicative and additive 
biases are consistent with zero at the sub-percent accuracy level . This agrees with the 
conclusions from fig.l and 2. Note that there are roughly 0.1% residual systematic errors, 
which are likely from third order corrections. On the other hand, the data on the tables show 
that the pixelation effect becomes less important for larger galaxies. This is consistent with 
our intuitions. In the next section, we include noise in the simulations to test our method 
under realistic conditions. 

4.3 Tests with Photon Noise 

There are three types of photon noise: 1. the spatial fluctuations of the sky background; 2. 
The Poisson photon counting noise of the background photons; 3. The source Poisson noise. 
The noise treatment defined in [25] only deals with the background noise. We currently do 
not have a way to treat the source Poisson noise, which can certainly introduce systematic 
errors to the shear recovery. In the weak lensing community, one typically use the signal- 
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Small Gal. 


Medium Gal. 


Large Gal. 


Largest Gal. 




= 1 


mi(10~ 3 ) : 2.0 ±1.5 


-0.8 ± 1.5 


3.8 ±1.5 


2.8 ± 1.4 






ci(10" 5 ) : 9.6 ±4.8 


4.3 ±4.8 


-2.1 ±4.7 


-4.8 ±4.7 


-Lpixel 


= 2 


3.2 ± 1.5 


-0.5 ± 1.5 


3.9 ±1.5 


2.9 ± 1.4 




9.4 ±4.8 


4.2 ±4.8 


-2.1 ±4.7 


-4.8 ±4.7 


Lpixel 


= 3 


7.2 ± 1.5 


0.6 ± 1.5 


4.2 ± 1.5 


3.0 ± 1.4 






10.2 ±4.8 


4.3 ±4.8 


-2.1 ±4.7 


-4.8 ±4.7 


-Lpixel 


= 4 


-2.5 ±1.5 


-1.9 ± 1.5 


3.5 ± 1.5 


2.8 ±1.4 




8.9 ±4.8 


4.1 ±4.8 


-2.1 ±4.7 


-4.8 ±4.7 


Lpixel 


= 5 


-91.2 ±1.3 


-30.6 ±1.4 


-4.6 ± 1.4 


0.4 ± 1.4 






9.8 ±4.3 


4.1 ±4.6 


-2.1 ±4.7 


-4.8 ±4.6 


Lpixel 


= 6 


-372 ± 1 


-174.5 ± 1.2 


-58.4 ± 1.4 


-16.2 ±1.4 




7.5 ±3.1 


3.2 ±3.9 


-2.0 ±4.4 


-4.7 ±4.6 



Table 1. The multiplicative bias n%\ and the additive bias c\ of the first component of the reduced 
shear measured using the method of this paper. In each data cell, the upper value is mi in unit of 
10 -3 , and the lower value is c\ in unit of 10 -5 . The results are shown for six choices of the pixel size 
and four choices of the galaxy size. Noise is not included in this set of simulations. 







Small Gal. 


Medium Gal. 


Large Gal. 


Largest Gal. 


Lpixel 


= 1 


m 2 (10~ 3 ) : -0.9 ±1.5 


-1.6 ± 1.5 


-0.4 ± 1.5 


1.6 ± 1.4 




c 2 (10- 5 ) : 7.3 ±4.8 


-0.5 ±4.8 


7.4 ±4.7 


1.2 ±4.7 


Lpixel 


= 2 


0.3 ± 1.5 


-1.3 ±1.5 


-0.3 ± 1.5 


1.7 ±1.4 




7.6 ±4.8 


-0.4 ±4.8 


7.5 ±4.7 


1.2 ±4.7 


Lpixel 


= 3 


4.3 ± 1.5 


-0.2 ± 1.5 


0.0 ±1.5 


1.8 ±1.4 






8.9 ±4.8 


-0.3 ±4.8 


7.5 ±4.7 


1.2 ±4.7 


Lpixel 


= 4 


-5.4 ± 1.5 


-2.7 ±1.5 


-0.7 ± 1.5 


1.5 ± 1.4 






10.2 ±4.8 


-0.3 ±4.8 


7.4 ±4.7 


1.2 ±4.7 


Lpixel 


= 5 


-100.3 ±1.3 


-32.2 ± 1.4 


-8.5 ± 1.4 


-0.5 ±1.4 




11.0 ±4.3 


-0.1 ±4.6 


7.4 ±4.7 


1.2 ±4.6 


Lpixel 


= 6 


-420.2 ±0.9 


-192.3 ±1.2 


-62.8 ±1.4 


-16.8 ±1.4 






54.2 ± 2.8 


4.4 ±3.9 


7.2 ±4.4 


1.2 ±4.6 



Table 2. Same as table 1, except that it is for multiplicative and additive biases (m^ and C2) of the 
second component of the reduced shear. 

to-noise-ratio (SNR) to denote the noise amplitude relative to the source signal within the 
halfdight radius of the source (post-seeing). Here, we point out that at least within the scope 
of our method, it is important to differentiate between background and source noise for the 
reasons just discussed. In other words, for any given SNR, whether most of the noise is due 
to the background or the source can significantly affect the shear recovery accuracy. The 
ratio of the background noise to the source noise is determined by the luminosities of the 
source and the background. 

The purpose of this section is to study the shear recovery accuracy under different SNR. 
For a given SNR, we consider three cases: 1. solely background noise; 2. the background 
and source Poisson noises are even; 3. solely source Poisson noise. The background noise 
in our simulations is generated as the Poisson noise of a homogeneous background. To save 
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Case 1 Case 2 Case 3 



Figure 3. Sample images of the same galaxy under 9 different noise conditions. Images in the same 
row have the same SNR, which is shown on the left. Images in each column share one of the three 
cases: 1. solely Poisson noise of the background; 2. even contributions from the Poisson noises of the 
background and the source; 3. solely source Poisson noise. 



computational time, the background fluctuations are not included in this paper. This should 
not affect our main conclusions. Fig. 3 shows sample images of a single simulated galaxy with 
three SNR values: 10, 20, 40. For each SNR, the three cases just discussed are shown in the 
same row. 

The galaxies and PSF used in the simulations of this section are discussed in §4.1. The 
pixel size L p i xe i is taken to be 4, which matches what is assumed in the GREAT08 project. 
For treating the background noise using the method of [25], we randomly generate an image of 
pure background noise for each source image. Similar to what is done in §4.2, when a galaxy 
is generated, it is used for different shear values and noise levels for the purpose of saving 
time. For each simulated galaxy, the background Poisson noise in each pixel is generated as 
a Gaussian random number with a unitary variance 4 , and rescaled according to the assumed 
SNR due to the background noise (SNR#). A separate image of pure background noise is 
generated in the same way with different random seeds. The source Poisson noise in each 
pixel is generated as a Gaussian random number with unitary variances multiplied by the 
square root of the source flux of the same pixel. It is then rescaled according to the value 
of SNR due to the source Poisson noise (SNR5). Since the source Poisson noise and the 
background noise are not correlated, we use the following equation to relate SNR, SNR#, 
SNR 5 : 

111 , , 

9 = 5- + 9 (4-3) 

SNR 2 SNR| SNR| V ; 

For each given SNR, case 1, 2, and 3 correspond to SNRb/SNRs = 0, 1, and 00 respectively. 
We consider six choices of SNR: 10, 20, 30, 40, 60, 80. The main results are shown in table 



4 The Gaussian distribution is a good approximation to the Poisson distribution when the background 
photon number in each pixel is not a very small number. 
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Case 1 


Case 2 


Case 3 


SNR= 


80 


mi(10- 3 ) : 0.2 ± 1.4 


5.0 ± 1.2 


8.8 ±0.8 






ci(10~ 5 ) : -7.1 ±4.7 


-6.9 ±3.8 


-4.6 ±2.6 


SNR= 


60 


0.7 ±1.8 


9.3 ±1.4 


16.6 ±0.9 






-8.4 ±5.9 


-8.2 ±4.6 


-5.1 ±2.8 


SNR= 


40 


1.9 ±2.6 


21.3 ±2.0 


39.2 ± 1.0 






-11.1 ±8.6 


-10.9 ±6.5 


-6.1 ±3.3 


SNR= 


30 


2.9 ±3.5 


38.0 ±2.7 


72.5 ± 1.2 






-13.9 ±11.4 


-13.8 ±8.6 


-7.2 ±3.9 


SNR= 


20 


5.0 ± 5.5 


87.5 ±4.2 


179.5 ± 1.8 






-19.9 ±17.8 


-20.5 ± 13.6 


-10.0 ±5.7 


SNR= 


10 


10.3 ±13.3 


451.9 ± 12.3 


1529.5 ±6.8 






-40.8 ±43.4 


-56.0 ±40.1 


-34.8 ±22.2 



Table 3. The multiplicative bias mi and the additive bias c\ of the first component of the reduced 
shear measured using the method of this paper. In each data cell, the upper value is m\ in unit of 
10 -3 , and the lower value is c\ in unit of 10 -5 . The results are shown for six choices of SNR. For 
each SNR. three cases are considered: 1. solely background noise; 2. the background noise and source 
Poisson noise have equal amplitudes; 3. solely source Poisson noise. 







Case 1 


Case 2 


Case 3 


SNR= 


80 


m 2 (W"'- i ) : -0.3 ± 1.4 


4.3 ±1.2 


8.3 ±0.8 






c 2 (10- 5 ) : 3.0 ±4.7 


2.2 ±3.8 


1.4 ±2.6 


SNR= 


60 


0.1 ± 1.8 


8.3 ±1.4 


15.9 ±0.9 






2.9 ±5.9 


1.9 ±4.6 


1.1 ± 2.8 


SNR= 


40 


0.8 ± 2.6 


19.6 ±2.0 


38.1 ± 1.0 






2.2 ±8.6 


1.2 ±6.5 


0.3 ±3.3 


SNR= 


30 


1.4 ±3.5 


35.7 ±2.7 


70.7 ± 1.2 






1.1 ± 11.4 


0.2 ±8.6 


-0.5 ±3.9 


SNR= 


20 


2.8 ±5.5 


84.0 ± 4.2 


176.5 ± 1.8 






-2.7 ± 17.8 


-2.7 ±13.6 


-2.2 ±5.7 


SNR= 


10 


6.7 ±13.3 


444.5 ± 12.3 


1516.7 ±6.8 






-26.2 ±43.4 


-23.7 ±40.1 


-15.3 ±22.2 



Table 4. Same as table 4, except that it is for multiplicative and additive biases (7712 and C2) of the 
second component of the reduced shear. 

3 and 4. 2.6 x 10 7 mock galaxies are used to measure the biases in every case of the tables. 

The results in tables 3 and 4 show that the shear recovery accuracy in our method 
strongly depends on whether the noise is mostly due to the background or the source. As 
shown in all the first cases in the tables, the shear measurement errors caused by pure 
background noise can be cleanly removed with the treatment of [25]. We can almost confirm 
the sub-percent level accuracy for SNR~ 10 in this case. The statistical errors, though 
systematically larger than those in other cases of the same SNR due to the nature of the 
noise treatment of [25] , can be further narrowed down to check the accuracy of the method 
at even smaller SNR with a larger galaxy ensemble. In the second and third cases of each 
SNR, the multiplicative biases are clearly all larger than their statistical errors, and reach 
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> 1% level for SNR< 60. To achieve a sub-percent level accuracy, we find that it requires 
SNRs > 80 — 100 in our method. It corresponds to collecting roughly 10 4 source photons per 
galaxy. Indeed, the limit on the SNRs is so far the only requirement for achieve very high 
accuracy in our shear measurement method. 

Physically, the shear recovery errors in case 2 and 3 are caused by the fact that the galaxy 
shapes are intrinsically modified by their own Poisson noise. It seems likely to smooth out 
the source Poisson noise if the half-light radius of the source covers enough pixels. However, 
it is hard to see how useful such a procedure is, because for a given SNRg, more pixels within 
the half-light radius mean a larger Poisson noise in each pixel. In the end, the total source 
flux should always has a relative uncertainty of order l/SNRg. Along this line of thinking, 
one can even guess a simple relation between SNRg and the amplitude of the multiplicative 
bias m through the following procedure: 



5g (S(f 2 )) 

9 (f 2 

(5(f)) ~ ((fs + In? - /I) ~ (2fsf N + fl) ~ (f 2 N ) 



m ~ — 77TT ( 4 - 4 ) 



m ~ 



</|> SNR| 



In the above estimation, g stands for the reduced shear, / refers to the total flux within 
the half-light radius of the source, and the subscripts S and iV denote the source and noise. 
Note that the (/jy) due to the background noise is statistically estimated using a neighboring 
map of pure noise and subtracted in our method. The remaining shear measurement errors 
therefore solely come from the source Poisson noise. The results in tables 3 and 4 indeed 
indicate that m ~ 6O/SNR5, except when SNR5 < 10. Note that this is true in both case 
2 and case 3. We caution that the derivation in eq.(4.4) is not rigorous. It only provides a 
possible understanding of the results seen in tables 3 and 4. We will study source Poisson 
noise more systematically in a future work. 



5 Summary 

Based on [24, 25, 38], we have established a robust way of measuring the cosmic shear to the 
second order in accuracy. The method is well defined regardless of the morphologies of the 
galaxies and the PSF. We have also provided a useful theorem for judging the accuracy of 
any shear measurement method at the second order based on its properties at the first order. 

For our method to achieve the accuracy at sub-percent level, the CCD pixel size is 
required to be not larger than about 1/3 of the FWHM of the PSF, regardless of whether the 
PSF has a power-law or exponential profile at large distances 5 . Using more than 10 mock 
galaxies of unrestricted morphologies, we have tested the accuracy of this method under 
different noise conditions. We find that it is useful to separately discuss the background and 
source noise for any given SNR. The background noise, which is uncorrelated with the source 
flux, can be removed in a simple and clean way using the method of [25] . In our simulations 
with only background noise, the shear measurement errors are found to be less than 1% 
for SNR as low as 10, and the conclusion can likely be extended to even smaller SNR with 
simulations of a larger galaxy ensemble. On the other hand, the source Poisson noise, which 

5 For PSFs with strong diffraction spikes, we need to further test the method. This will be done in a future 
work. 
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strongly couples with the distribution of the source flux, remains to be the main cause of the 
shear measurement errors in our method. For a sub-percent level accuracy, we require the 
SNR of the source Poisson noise to be > 80 — 100. This corresponds to collecting about 10 4 
source photons per galaxy. The treatment of source Poisson noise is unclear at present, and 
will hopefully be addressed in a future work. 
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Appendix — The Relation Between Two Definitions of Cosmic Shears 

The convention of defining the cosmic shear /convergence in this paper is different from what 
is used in [24, 25, 38]. More specifically, as shown in eq.(2.1), we define shear /convergence 
using the following formula: 

( xf \ = ( 1 - k - 7i -72 
V x 2 ) V ~72 1 - k + 71 

In [24, 25, 38], however, shear /convergence is defined as: 

(x% \ = ( 1 + k' + i x 72 
\x\) V 7 2 l + «'-Ti 

Note that these two conventions are equivalent up to the first order in shear/convergence, 
but not to the second order. In other words, the values of 71, 72, and n are equal to those of 
7i, 72, and k' respectively only when the second and higher order lensing terms are neglected. 
It is straightforward to derive the relation between the two conventions based on the following 
identity: 

/ 1 + k' + Vj 72 
^ 72 1 + «' - tJ 

To the second order in shear/convergence, we get: 

/ ,2,2,2 
K = K+ K +7l+72> 

7 1 = 71(1 + 2k), 

7 2 = 72(1 + 2k). (5.4) 




\_/l-K-7i -72 \ 
/ I -72 1 - « + 7i / 
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As a result, we find: 

7i )2 (l - «') = 71,2(1 + K ). (5.5) 

Therefore, for theorists intending to make second order predictions for the cosmic shears, con- 
ventions should be explicitly mentioned, as they carry different meanings and consequences 
at the second order level. 
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